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We represent QCD at the hadronic scale by means of an effective Hamiltonian, H , 

formulated in the Coulomb gauge. As in the Nambu-Jona-Lasinio model, chiral 

Q , i symmetry is dynamically broken, however our approach is renormalizable and also 

fl} ■ includes confinement through a linear potential with slope specified by lattice gauge 

rf\ ' theory. We perform a comparative study of alternative many-body techniques for 

approximately diagonalizing H: BCS for the vacuum ground state; TDA and RPA 
l/"") t for the excited hadron states. We adequately describe the experimental meson 

and lattice glueball spectra and perform the first relativistic, three quasiparticle 

calculation for hybrid mesons. In general agreement with alternative theoretical 

approaches, we predict the lightest hybrid states near but above 2 GeV, indicating 

^. ' the two recently observed J PC = 1 I" exotics at 1.4 and 1.6 GeV are of a different, 

perhaps four quark, structure. We also detail a new isospin dependent interaction 
Qv , from qq color octet annihilation (analogous to ortho positronium) which splits I = 

I ■ and 1 = 1 states. 

o\ 

(—} 1 Introduction and Model Hamiltonian in the Coulomb Gauge 

In the past several years we have implemented Iffl'Bcl an ambitious QCD pro- 
i-G ' gram to comprehensively investigate hadron structure. Central to this project 

is a realistic effective Hamiltonian which we develop from QCD through renor- 
O !* malization and subsequently diagonalize using established many-body tech- 

niques. We are particularly interested in the signature QCD prediction of ex- 
otic states (hadrons with quantum numbers not possible in pure qq or qqq mod- 
els). Early evidence a for the existence of an exotic hadron with J PC = l r T + , 
denoted p, has finally been confirmed with the recent BNL observation El of 
two states at 1.4 and 1.6 GeV. These measurements, along with anticipated 
future hybrid experiments (COMPASS at CERN and Hall D at JLab), have 
motivated the study detailed in this paper. 

The theoretical situation is not satisfactory: lattice results U indicate that 
the lightest reliable hybrid meson mass with quantum numbers 1 *~ should 
appear at (2.0 ± 0.2) MeV. On the other hand, vintage bag model results □ 
and a recent QCD sum rule calculation El agree in predicting exotic masses, 
between 1.4 and 1.7 GeV. The sum rule predictioncl claims a rigorous 1.9 GeV 
upper bound which is clearly below the lattice result. The lattice calculations, 
however, do require extrapolating for the u, d sector and are therefore less 
accurate for the light quark hadrons. Hence, further insight is definitely needed 
to determine if these exotic states are predominantly qqg objects or belong to 
other Fock space sectors, such as qqqq. 



Our many-bocbz Hamiltonian approach, which has been successfully tested 
for both glueballsEl and conventional mesons □, addresses this issue. Before 
treating the 3-body hybrid meson spectrum we present our model Hamilto- 
nian in detail and briefly summarize our previous results. Xiur QCD inspired, 
effective Hamiltonian is formulated in the Coulomb gauge E£l 



H= /"<ix* t (x)(-m-V + /3m)*(x)+rr fdx{U a -U a + B^-B^) 

-i|dxdyp a (xM|x-y|K(y), (1) 



containing both quark, \I>, and gluon, A a , B^ = V x A a , fields with color den- 
sity p a — ^t^j-* + f abc A b ■ II C . We adopt the standard current quark mass, 
to, values for the u, d, s and c flavors; m u = rrid = 5 MeV, m s = 150 MeV 
and to c = 1200 MeV. The linear interaction, Vl = err, is obtained from lat- 
tice measurements and Regge phenomenology yielding a — 0.18 GeV 2 . For 
certain observables we supplement this with the canonical Coulomb potential 
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Vc = — ^ (ct s = f* - ~ 0.2 — 0.4). The normal mode field expansions are 



*(*) = E 



A ^=Ji£ 



n?(x) 
rfk 



(2tt) 3 



dk 1 



„ik-x 



=[o?(k)+o? T (-k)]e* 



(^3yf[<(k)-«f(-k)]e^ 



w cA (k)6 cA (k)+u cA (-k)4(-k) 



ik-x 



(2) 



and with the Coulomb gauge transversality condition, k-a a (k) = 0, yields the 
commutation relation 



K(k),of (q)] = 5 ab (2n) 3 6^(k- q)(% - kk) ■ 



(3) 



We note the color density-density interaction,, n a (x) V^ ( | x — y | ) p a (y ) , only pro- 
duces color singlet states in the spectrum bOI3 when V is a linear, harmonic 
oscillator or any Fourier transformed confining interaction with a singular be- 
havior as q — ► k. This is a stringent check in our calculations, since the 
matrix elements must always contain factors of q — k to cancel this singular- 
ity. In particular, for our chosen linear potential, we have in momentum space 
V-(|q-k|) = -^. 



2 Ground State and BCS Mass Gap Equations 

We now wish to solve H^f = E^ as accurately as possible. In this section we 
focus on the ground state and introduce the Bardeen, Cooper and Schrieffer 
(BCS) transformation. We proceed by normal-ordering the Hamiltonian using 
the basis, Eq. (Eh, and minimize the ground state, or vacuum, expectation 
value. The key concept in this approach is that of a quasiparticle. The bare, 
current operators a, b, d are rotated to improved quasiparticle operators a, B, 
D by means of the BCS transformation. The vacuum state, |fi), is determined 
by the relations B\fl) = D\fl) = a\il) = 0. The gluon and quark BCS rotations 
are 



a?(k) = coshe fe <(k) + smh0 fe af (-k) (4) 

B c \(k) = cos— -6 cA (k) - Asin— 4 A (-k) 

D c \(-k) =cos— d c \(-k) + A sin— b f cX (k) , 

where Q k , Ok/2 are the BCS angles, further specified below. This also intro- 
duces a redefinition of the spectral function, Wk, and counter-rotated quasi- 
particle spinors U , V such that the expansions of the fields remain formally 
invariant. 

The specific variational parameters are the quark gap angle, <pk, related 
to the BCS angle by tan((j>k — &k) — m/k, and the gluon self-energy, Wfe, 
satisfying io k = ke~ 2Bk . The vacuum contains correlated Cooper pairs and 
explicitly breaks chiral symmetry due to the form of u> and <f> bOEj. The 
functional variation with respect to the two parameters 6 k , Q k 

generates two integral (mass gap) equations for the quark and gluon sectors, 
respectively 

k sin cj) k -m cos 4> k = - / V( \ k - q| ) [sin <fa cos <§> q k-q - sin <\> q cos <j> k \ (5) 

6 J (Zir)* 

o _ - _ 3 f dq y ... ,~ , 2 . (w 2 q -w 2 ^ 



C^-^^Gk-qlXl-Kk.qn l^-JL . (6) 



Significantly, the BCS vacuum is stable against quasiparticle pair creation since 
the "anomalous" terms in the Hamiltonian of the type aa, a^a^ , BD, B^D^ 
are also cancelled by the same rotation. 

For the linear potential 1/q 4 , simple dimensional analysis reveals Eq. (ra) 
is UV finite while Eq. (0) is logarithmically divergent. Hence we impose a 
momentum cutoff A = 4 GeV and proceed to solve both numerically after per- 
forming a standard three dimensional reduction. The details of this calculation 
are given elsewhere u. We verify that our variational solution has a minimum 
in energy (traditional Mexican hat shape). 

3 TDA and RPA for Mesons and Glueballs 

With these quasiparticle degrees of freedom we now construct the excited states 
from this vacuum. For both the quark and gluon sector we represent hadrons 
as quasiparticle pairs (eg. qq for mesons and gg for glueballs) and angular 
momenta couple to form states of good J PC . Next we invoke the TDA at the 
lp-lh level and diagonalize the model Hamiltonian in this truncated space. 
A subset of our extensive TDA calculations is displayed in Fig. [l] for the 
pseudoscalar mesons. In general there is broad agreement with the data, the 
most notable exception being for the 7r and rj. 

The insufficient w 200 MeV, it/ p mass splitting and related issue of chiral 
symmetry motivated our improved RPA treatment. Now the pion creation 
operator generalizes to 



)t = yy v „v 



x mW 3 - Y ijQiQ : 



with pion state \n) = Q'\RPA) and improved vacuum satisfying Q\RPA) = 0. 
Here qi, q t are again the BCS rotated quasiparticle operators and Xij, Y^ 
are the RPA wavefunction components obtained from the coupled equations 
of motions generated by 

(tt\[H,Q^\RPA) = M n (n\Q i \RPA) . (7) 

In the chiral limit (to = 0) the chiral condensate operator commutes with 
Q' and we rigorously compute M v = consistent with Goldstone's theorem. 
For to = 5 MeV explicit chiral symmetry breaking yields M^ — 294 MeV, 
significantly better than the TDA value (note, the physical pion mass is repro- 
duced with to = 3 MeV) . Since both TDA and RPA produce comparable spin 
splittings, we therefore conclude that chiral symmetry, which only the RPA 
respects, is responsible for most of the ir/p mass difference. 
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Figure 1: TDA (dots), RPA (dashes) and data (bars) for the light pseudoscalar and vector 
meson spectra. 



Similarly, we generate the glueball spectrum for two quasiparticle gluons. 
Now chiral symmetry is not an issue and the RPA and TDA spectra agree to 
within a few percent. Our calculations, using the same string tension as above, 
reproduce the lattice measurements as discussed in Confinement II. 

A fundamental test for QCD is the existence of exotic mesons (quantum 
numbers not possible in artv qq model). In particular, it has been speculated 
that two recently observed^ states with isospin 1 and J PC = 1 v at 1.4 and 
1.6 GeV contain explicit glue. Of several gluonic scenarios, glueballs (oddballs) 
must be eliminated since they have isospin 0. We therefore investigate in the 
next section if these 1 ^ states are hybrids. 



4 Three-Body TDA for Hybrid Mesons 

We formulate the hybrid meson as a 3-body problem (qqg) with TDA color 
wavefunction 



\hybrid) = BW"a ! \BCS) = [[B 



it 



■2^1 



a%\BCS), 



where the three BCS quasiparticles have momenta q, q, and g referred to 
the hybrid cm frame. The optimal choice of relative variables is q+ = ^rp , 



q_ = q— q which facilitates implementing rotational, parity (P) and C-parity 
symmetries. The complete TDA wavefunction ansatz with spin (A) and color 
(c, a) indices is 

\hybrid)= £ //^^^^(q + ,q-)^< c (q)^(q) a f( g )|^). 

\ q \-Xcca J V ' V ' 

This wavefunction's rotation properties are governed by the three tensor 
indices X q , Xg, A, and the two arguments q+, q . To construct an appropriate 
SU(2) representation we need to couple five angular momenta (three quasipar- 
ticle intrinsic spins and two orbital, L±, associated with q±) to give a total 
J, raj. To accomodate the Coulomb gauge transversality condition we utilize 
a modified LS coupling scheme involving two intermediate angular momenta 
1 = L_|_ + 1 and L = 1 + L . The resulting angular momentum decomposition 
is 

F/^(q + ,q_) = ]T F^XlsQi+I k-l) ^(q+J^/G-) 

IL — L+LSm+m— 

(— l) A *(£ + m+l - \g\lmi) {L_m-lmi\LmLl 
; 2 Ag 2 



<o\o " >>q\Sm s )(-l)i- x «(Lm L Sms\Jmj} 



Due to the transversality condition mentioned above, q+-« = 0, so 

(im+iA|oo)yr + (q+)4 = o (8) 

and therefore we cannot couple the angular momenta of the gluon to give the 
intermediate state \lmi) = |00). 

With this coupling scheme the total parity is a product of the intrinsic 
parities of the quark- antiquark pair (-1), the intrinsic parity of the gluon (-1), 
and the two spherical harmonics yielding 

P = (-l)-(-l)-(-l) L + +L - = (-l) L + +L - . 

Similarly the total charge conjugation is a product of the quark-antiquark 
pair ((— \) L ~ + from equivalence to the exchange of all qq quantum numbers) 
and the gluon (-1), which is its own antiparticle and has odd C-parity. Hence 



C=(-l)-(-l) L - +s = (-l) 



1+L-+S 



Notice the extra C-parity sign which now permits exotic quantum numbers. 
We are only interested in the lightest hybrids which will have S and P waves 



Table 1: Possible Quantum Numbers for a qqg Hybrid Meson 
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Figure 2: Light and charmed exotic 1 + hybrid mesons. The TDA calculation is comparable 
to the lattice results but both disagree with observed resonances. 



for either of the angular momenta L + or L-. The resulting possible quantum 
numbers are displayed in Table nj. 

In principle all the states with the same J quantum numbers can mix 
but we do not address this issue here. We also note that our wavefunction 
normalization is not standard due to the transversality relation, Eq. (0) , which 
introduces an extra angular momentum coefficient. 

Finally, in analogy with our two-body TDA treatment the hybrid equation 
of motion is 

(qqg\H\hybrict) = E(qqg\hybrid) . (9) 

Using Eqs. ([D, 0) and the hybrid wavefunction expansion, we obtain the TDA 
equation of motion which now is a non-local equation in the two variables 
q + , q (12-dimensional problem) precluding practical matrix diagonalization. 
Therefore we utilize alternative ansatze for the radial wavefunction and varia- 
tionally evaluate the Hamiltonian matrix element. 

Separating the cm variables, the integrals can be reduced to 9-dimensions 
which we evaluated by Monte Carlo methods utilizing Lepage's code VEGAS. 
We employed an increasing number of points, up to several million, until sat- 



isfactory convergence was achieved. The angular wavefunctions are explicitly 
coded in terms of spherical harmonics with all magnetic spin sums and angular 
integrals evaluated numerically. The integrable IR singularity from our poten- 
tial received special attention. We divided the complete integral into different 
parts depending upon the Hamiltonian components and then used a change of 
variables, always placing the singularity at the origin, with a Jacobian trans- 
formation to concentrate points there. This is detailed in our previous work. 

In Fig. |2| we present the main result of this work, our prediction of the 
light and charmed exotic hybrid mass. We also compare to alternative lattice, 
flux tube and even vintage bag model approaches as well as the recent E852 
Brookhaven data. Notice that, with the exception of the dated bag model 
results, the different theoretical models generally agree but disagree with ob- 
servation. 

Because of the potential ramifications of this disagreement between mea- 
surement and predictions we further examined and tested our model and cal- 
culation in detail. First, to ensure our variational method is accurate, we 
compared several different radial trial wavefunctions. In particular, we used 
gaussian and also more complicated functions generated by the exact numerical 
solutions to simpler, 2-body problems (p meson, J/ <I>, glueball) with various 
potentials. We also added nodes as needed. We found the maximum variance 
(sensitivity) to wavefunction choice to be about 50 MeV in the hybrid mass. 
Second, we also considered whether chiral symmetry affects the hybrid meson 
spectrum and variationally formulated the RPA as above, but now with an 
additional constituent TDA gluon. Interestingly, we find RPA and TDA to be 
numerically equivalent. The reason is that the quarks are now in a color octet 
state but the chiral octet charge 



Ql = Jdx^(x) l5 T a *(x) 



is not conserved since the commutator with the Hamiltonian does not vanish. 
Further, there is also no Gell-Mann-Oakes-Renner relation and no Goldstone 
hybrid boson in the chiral limit m — > 0. Therefore we do not expect any major 
numerical effect in the hybrid spectrum, analogous to the 400 MeV downward 
shift of the chiral pion. This was confirmed by an explicit numerical calculation 
in a selected channel. 

Consequently, we are reasonably confident in our prediction that the light- 
est exotic hybrid mass is above 2 GeV and concur with other contemporary 
theoretical model results. This also strongly indicates that the observed exotic 
1 h at 1.4 and 1.6 GeV are not hybrid mesons but have an alternative, four 
quark or meson molecular structure. 
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Figure 3: TDA theory for conventional (cc) and non-exotic (ccg) states compared to the 
observed 1 J/ip spectrum 1 4. 



Our last key result concerns the charmonium J/^S spectrum and is shown 
in Fig. 3. By simply calculating D wave TDA cc states we have resolved the 
historic "overpopulation" problem. Ironically, if we then include our predicted 
non-exotic charmonium states, also shown in the figure, we now have an "un- 
derpopulation" problem which is more typical of confrontation between the 
conventional quark model and experiment. It would now appear that simply 
counting states may not be that effective in identifying gluon rich states. 

Finally, since any three quasiparticlc problem can be treated variationally, 
we have explicitly constructed valence baryons (qqq) states and evaluated both 
the nucleon (TV) and A 3 / 2 with our same Hamiltonian. Our A mass is in 
agreement with data, but only about 1/3 of the vV-A mass splitting is obtained. 
The other interesting case is the triple gluon glueball, ggg, whose study is still 
preliminary and will be reported in the near future. This exhausts the possible 
three body states which can form color singlets. 
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5 Conclusions 

In summary, the BCS, TDA and RPA many-body treatments are power- 
ful, effective methods for investigating hadron structure. Using a relativis- 
tic, field theoretical Hamiltonian with standard quark masses and only one 
(pre-determined) interaction parameter, our approximate many-body solutions 
yield reasonable descriptions of the vacuum, mesons and glueballs. In particu- 
lar, the RPA properly incorporates chiral symmetry and provides new insight 
into the condensate structure of the vacuum and chiral governance of the pion. 
Most significantly, our hybrid meson mass prediction is above 2 GeV and in 
reasonable agreement with both lattice and flux tube results. This strongly 
suggests that the recent observed exotic states below 2 GeV are not hybrids 
but more likely 4 quark states. Because our method is directly amenable 
to including higher Fock state components, we are currently calculating the 
molecular (qqqq) meson spectrum. Future work will also entail extensions to 
the nucleon strangeness content (pentaquark systems) and even dibaryons (six 
quark systems) which are both tractable in our approach. 
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